Distribution of velocities in an avalanche 
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For a driven elastic object near depinning, we derive from first principles the distribution of 
instantaneous velocities in an avalanche. We prove that above the upper critical dimension, d > d uc , 
the n-times distribution of the center-of-mass velocity is equivalent to the prediction from the ABBM 
stochastic equation. Our method allows to compute space and time dependence from an instanton 
equation. We extend the calculation beyond mean field, to lowest order in e = d uc — d. 

PACS numbers: 68.35.Rh 



Obtaining a quantitative description of the dynamics 
during an avalanche is of great importance for systems 
whose dynamics is governed by jumps, such as magnets, 
superconductors, earthquakes, the contact line of fluids, 
or fracture [THS]. In particular the motion of domain 
walls (DW) in magnets is important for many applica- 
tions, such as magnetic recording. It can be measured 
from the Barkhausen (magnetization) noise [SJ [7] , which 
is a complicated time-dependent signal. Its origin is due 
to an interplay between quenched impurities and the elas- 
tic deformation energy which tend to pin the DW, as well 
as the driving and magnetostatic forces. 

A major step forward was accomplished by Alessandro, 
Beatrice, Bertotti and Montorsi (ABBM) [§. who intro- 
duced, on a phenomenological basis, a stochastic equa- 
tion approximating the DW motion by a single degree of 
freedom. Although a crude description, this model has 
been used extensively to compare with experiments on 
magnets, with success in some "mean-field like" cases, 
and failure in other [9]. However, no microscopic foun- 
dation for the validity of this model exists. 

On the other hand, sophisticated field theoretic meth- 
ods were developed in the last decades to study systems 
with quenched disorder. In particular, for elastic inter- 
faces, relevant to describe DW motion, functional RG 
methods (FRG) [U H0lfl2] have recently allowed to derive 
the distribution of quasi-static avalanche sizes [Hfl HI] . 
Until now however, no description of the dynamics dur- 
ing an avalanche was available. In fact, since it involves 
much faster motion than the average driving velocity, it 
led to difficulties in the early FRG approaches [TT] . 

The aim of this Letter is to show how to compute from 
first principles the distribution of instantaneous velocities 
in an avalanche. We study a single elastic interface, of 
internal dimension d (total space dimension is D = d+1) 
at zero temperature, near the depinning threshold. The 
method works in an expansion around the upper criti- 
cal dimension d uc , with d uc = 4 for standard elasticity, 
and d uc = 2 in presence of long-range elasticity, e.g. aris- 
ing from dipolar forces. Remarkably, we find that for 
d = due (and above) and in the scaling limit, the n-time 
probability distribution (with n arbitrary) of the center 
of mass of the interface is equivalent to that of the ABBM 
stochastic equation, in terms of renormalized parameters 



which in some cases can be estimated. The two meth- 
ods are rather different in spirit, and the identification 
non-trivial. Our result establishes the universality of the 
ABBM model for d > d uc . In addition it allows to resolve 
the spatial structure, and gives the corrections to ABBM 
for d < d uc . 

Here we sketch a very simple derivation, for details and 
various subtleties involved we refer to |15j . Consider the 
equation of motion, in the comoving frame, for the local 
velocity of an interface driven at velocity v: 



(r]od t - Vl)u xt = d t F(vt + u xt ,x) 



m 2 u xt 



(1) 



It is obtained by time derivation (noted indifferently u 
or d t u) of the standard overdamped equation of motion. 
Here x is the d-dimensional internal coordinate, vt + u xt 
the space and time dependent displacement field and rjo 
the friction. F(u, x) is the quenched random pinning 
force from the impurities, with e.g. Gaussian distribution 
and variance F(u,x)F(u',x / ) — S d (x — x')A (u — u'). 
m 2 is the strength of the restoring force —m 2 {u x t — vt) 
(i.e. the mass, or spring constant), which flattens the 
interface beyond a scale L m ~ 1/m. In the small m, large 
L m , limit, studied here, the interface has the roughness 
exponent £ of the depinning transition, with u ~ x^ for 
x < L m and u ~ for L > L m . For simplicity we chose 
standard elasticity ~ V 2 , but it can be replaced by an 
arbitrary elastic kernel as needed in applications [U EH E] • 
Near the depinning transition, i.e. at small v, the inter- 
face proceeds via avalanches. This is easiest seen in the 
center-of-mass position u t = L~ d J u xt . There is a well- 
defined quasi-static limit v = + where u t — u(w), with 
w = vt the well position. The process u(w) jumps at dis- 
crete locations Wi, i.e. u(w) — L~ d J2i Si8(w — Wi), with 
Si the avalanche sizes. Their statistics was predicted via 
FRG, and checked numerically [RJl QH [17] . There, the 
bare disorder correlator A (it) flows, under coarse grain- 
ing, to the renormalized one A (it), which, at the depin- 
ning transition exhibits a linear cusp — A'(0 + ) > 0. This 
cusp is directly related to the moments of the normalized 
size distribution P(S), via [T4] 
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FIG. 1: Schematic plot of the instantaneous velocity (divided 
by v) as a function of vt for different v. The area under the 
curve is the avalanche size hence is constant as v — > + . The 
quasi-static avalanche positions Wi are indicated. 



S m ~ m~( d+< ») is the large-scale cutoff of P(S). Here we 
study the dynamics inside these avalanches, which occur 
for small v on a time scale r m ~ L z m <C Aw/v, where 
Aw is the typical separation of avalanches in the same 
space region, and z the dynamical exponent. Hence we 
are considering small enough v so that avalanches remain 
well separated, a condition equivalent to L m <C where 
is the standard critical correlation length jTUJ HI] near 
depinning (for m = 0) . This is illustrated on figure [T] 

The information about the dynamics in an avalanche 
is contained in the n-times cumulants C n — ut 1 ■ ■ ■ iit n , 
n > 2 (with u t = 0). In the limit v — >• + the prod- 
uct ii tl . . . ut n vanishes unless all times are inside an 
avalanche. The probability that exactly one avalanche 
occurs in a time interval T < Aw/v is p^vT 1 with 
po = L d /(S) the avalanche density per unit w. C n 
is thus 0(v), rather than 0(v n ), the hallmark of a 
non-smooth motion. In addition, C n obeys the sum 
rule L nd Jj_ 



-T/2,T/2]' 



dti... dt n u tl . . . u tn = p vT(S n ) 



0(v 2 ). It can be computed perturbatively in the (renor- 
malized) disorder. For n — 2 and to lowest order one 
finds 



u tl u t2 



L- d A'((T 
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where here and below n is the renormalized friction |16j . 
Integrating over time, one recovers ([2]). 

To obtain all moments at once, as well as the velocity 
distribution, we now compute the generating function 



Z[X] = L- d d v e^t x ^+^t) 



u=0+ 



(4) 



tained from the dynamical action S = Sq + Sdis of |T]): 
~ 1 m 2 )u xt (5) 



5*0 = / u xt (r)d t - V 

I xt 
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u x tu x t'd t dt'A(v(t - t') + u xt - u xV ) (6) 
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This yields 

Z[X] =L- d d v J 2?[u]2?[u] e - s +/x t ^(^+^) 

with Z[0] = 0. We write 

d t d t 'A(v(t - t') + u xt - Uxt') 
= (v + u xt )d t < A' (v(t - t') + u xt - u xt ') 
= (v + M xt )A'(0+)c^sgn(t -*') + ... 



«=0+ 



(7) 
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where we have used that the interface is only moving 
forward (Middleton theorem [IS]). We can thus rewrite 

. ., where 



the disorder term as S = S^° c 
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is the so-called tree- level or mean-field action [16] . The 
terms neglected are O(A"(0 + )) and higher derivatives, 
and we have shown that they contribute only to O(e) to 
Z[X], hence can be neglected at tree level. 

We now study the tree approximation for Z[X], i.e. ^ 
with S'dis replaced by ^ . Thus the highly non-linear ac- 
tion ^ has been reduced to a much simpler cubic theory! 
Even more remarkably, u x t appears only linearly in (|9|, 
and viewing u as a response field, the tree level theory is 
equivalent to the following non-linear equation: 



m 2 )u xt - A'(0+)u 2 xt + X xt = (10) 



We denote u xt the solution of this equation for a given 
source X x t- Performing the derivative w.r.t v in ^ gives 

Z[X]=L- d ( X xt -A'(0+)(^ xt ) 2 (11) 

Jxt 



L- d [ (- V d t -W 2 x +m 2 )u x xt =m 2 L- d f 



u x 



where we have used equation ( 10 1 and, in the last equal- 
ity, assumed that u x vanishes at large t and x. To an- 
alyze the result, it is convenient to use dimensionless 
equations, replacing x — > x/m, L — > L/m, t — > T m t, 
v -> vv m , X -> X/S m and ii xt -> u xt /m 2 S m , where 
v m = S m m /r m , and r m = r\jm 2 . From now on we 
use these units, and consider the center-of-mass velocity, 
thus choosing X xt = X t uniform. 

The 1-time probability at time t = is given by At = 
XS(t) through its Laplace transform 



The average over disorder (and initial conditions) is ob- 



Z(X) = L- d d v e Ldx ^+u) 
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it = u t =o and the notation Z reminds us that we use 
dimensionless units. u xt = u t and we need to solve 



(d t - l)u t + u 



r.2 



with u t — > at t — ±00: 



ut 
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Inserting into (12 1 gives 



Z{\) 



ln(l - A) 



(13) 



(14) 



(15) 



Calling Tj the duration of the z-th avalanche out of N, and 
defining (r) := ^ J2< the mean duration, the probabil- 
ity Pa that t = belongs to an avalanche is p a = pqv(t). 
Hence the total 1-time velocity probability is P(u) — 
(1 — p a )5(v + u) + p a P(u) where P(u) is the probability 
given that t = belongs to an avalanche. Both P and P 
are normalized to unity. One notes the two (always) ex- 
act relations (u)p — 0, p a (u + v) p = v. Hence for v — + 
one has pq(t)(u) p — 1 and, in dimensionfull units Z(X) = 



-Z(m d v m \) = L- d p Q ( T ) J d«P(u)(, 



1). We 



thus obtain, in the slow driving limit, the instantaneous 
velocity distribution in the range v <C u ~ w m («o being 
a small velocity cutoff): 



PoiVVf,, \v m 



p(x) = —e 
x 



We defined S m = (m£) d w m = L d S m /T m . Hence 
(«) p ~ w m /ln(^-). Note that (i) p(x) is not a prob- 
ability, but is normalized by J dxxp(x) — 1 (ii) the 
quantity which is distributed according to p{x) is x = 
T m f Uxt/S m , which does not contain the factor L~ d . 
Similarly one obtains the n-time distribution of 



the center-of-mass velocity solving ( 13 ) with A t = 
ELi V(* - *,0. noting Zlj := 1 - e -l*=W T » 



z„(A 1 ,...,A„)=-in( ]r ni-\} n 2 

\AC{1 n}i£A {i,j}CA,i<j 

(17) 

For n = 2 one finds Z 2 = — ln(l — Ai — A2 + AiA 2 z) with 
2 = 1 — e"'* 2 "' 1 '/ 1 "™. From this we obtain (i) the proba- 
bility qi2 = vq'12 that both t\ and <2 belong to the same 
avalanche and the velocity distribution P conditioned to 
this event: 



p(ui,u 2 ) = ^ _ n . hi — 1 — 

(1 - e ^y^Iui V 1 - e / 



(18) 



(19) 



with t = — ix|/T m , q'^m = ln(l/z), and Ji(ar) is the 
Bessel-7 function of the first kind. The probability that 
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FIG. 2: "Pulse-shape": The normalized velocity at time t 
in an avalanche of duration r for r <C T m (lower curve) to 
t 3> T m (upper curve). 



ii but not *2 belongs to an avalanche is 
, - . 1 / Ui \ 
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(20) 



with p' Q = q[ + q' 12 . Since the probability that there 
exists an avalanche starting in [t\, t\ +dti] and ending in 
[h, * 2 +rf*2] is —dtidt2dt 1 dt 2 qi2 we obtain the distribution 
of durations t as 



P(r) = 



1 



,—t/t„ 



P0V„ 



(1 



r/r m )2 



(21) 



(16) For small durations r <C r m , P(r 
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ment with the above, using ln( — ) « ln( 



— J 1 — a, cut off at 

PQVmT 

In (7^) m good agree- 
Note that 



r / \ V 

q' 12 P(0 + ,0 + ) is proportional to the probability that an 
avalanche starts at t\ and ends at £2- 

The "shape" of an avalanche with duration r can then 
be extracted from the probabilities at 3 times (ti, i 2 , £3) = 
(0, t, t) setting ii\ = ii 3 = + . From the generating func- 



tion (17 1 for 3 times, the probability distribution for the 
intermediate-time velocity is P{u2) = b 2 U2e~ U2b , with 
v m b '■= — — h y~ — 1 resulting in the average "shape" 



U 2 = r = V r , 

b 



4smh (24 



sinhf^) 



(22) 



This interpolates from a parabola for small t < r m to a 
flat shape for the longest avalanches (see Fig 2.). This 
result holds for an interface at or above its upper critical 
dimension, which previously was used [7] on the basis of 
the ABBM model. 

We now clarify the relation to the phenomenologi- 
cal ABBM theory [5]. The latter models the inter- 
face as a single point driven in a long-range correlated 
random-force landscape, F(u), with Brownian statistics. 
It amounts to suppressing the space dependence in ([!]), 
hence corresponds in our general model to the special 
case d = and Ao(0) — Ao(u) = a\u\. The instanta- 
neous velocity v = u t + v satisfies the stochastic equation 



4 



■ndv = m 2 (v — v)dt + dF where dF 2 
ciated Fokker-Planck equation 



2ervdt , with asso- 



rt Q = d v 



-d v (\jQ) + m 2 (v-v)Q 



V 



(23) 



for the velocity probability Q = Q(v, tjvj, 0). For 
v > it evolves to the stationary distribution Qo( v ) — 
v - v /^ v v/v m -t e - v /v m /r(v/v m ) with v m = S m /r m and 
here S m = a/m and T m = rj/m 2 . For v = + one recov- 



ers ( 16 ), up to a normalization which entails a small-scale 

+ one finds the propagator 
-1-^,0) with 



cutoff. Similarly for v 
g(v,<|v 1 ,0)=i;- 1 Q(^ 



Q(v 2 ,t\v!,0) =v ie v 



s 1 — 
p{vi,v 2 )-\ e 

Vl 



(24) 



and p{yi, V2) given in Eq. ( 19 ). Q{v%, t\v\, 0) is solution of 
(23) with Q(v 2 ,0 + |vi,0) = 5(v 2 -vi). The piece - 5(« 2 ) 
corresponds to avalanches which have already terminated 
at time t, and is necessary for Q to conserve probabil- 
ity. The joint distributionQ(v2,t\vi,0)^-e~ Vl reproduces 



the 1-time and 2-times probabilities given in Eqs. ( 18 ) 
and (20), up to a global normalization. More gener- 
ally, since v(t) is a Markov-process, the n-time veloc- 
ity probability obtained from (10 1 is q' lp P(ui, . . . ,u n ) = 
iliijtj). 



Several remarks are in order. The first one is specific 
to the ABBM model: Since it is the zero-dimensional 
limit of ([!]), the dynamical-action method can be ap- 
plied. Hence we just found that for the ABBM model 
at v — + the tree approximation is exact. In the field 
theory it means that the effective action T equals the 
bare action S, and there are no loop corrections. Hence 
A'(u) = Aq(u) = — ersgn(w) is an exact FRG fixed point 
(with £ = 4 — d) as noted in [T3] . Crucial for this remark- 
able property is that the force landscape is a Brownian, 
and even in d = 0, this is not valid for any other, e.g. 
shorter ranged, force landscape. In that sense, the model 
proposed by ABBM [5], although unnatural from a mi- 
croscopic point of view, appears extremely judicious. 

Second, since a realistic interface in a short-ranged ran- 
dom force is described for d > d uc by the tree approxi- 
mation, we proved that the temporal correlations of its 
center-of-mass velocity for v — > are given by the ABBM 
model. Only two parameters enter, rj and S m , which in 
d = 4 acquire a logarithmic dependence on m |14j . 

Third, it is not expected that S = T extends to finite 
driving velocity v > 0; hence whether the phenomenol- 
ogy of the ABBM model with an avalanche exponent r 
dependent on v has anything to do with realistic interface 
motion remains an open question. 

Fourth, the present theory allows to go beyond the 
ABBM model in several ways: In d > 4, the non-linear 
equation ( 10 1 allows to study the full time- and space- 



dependence of velocity correlations, as was done for the 
statics in [T3]. Second, including loop-corrections al- 
lows to compute corrections in a systematic expansion 



in d = 4 — e [T5] . The main result for the 1-time velocity 
distribution for vq <C v <C v m is to first order in e 



P(v) - l/v a 



l-e(l-Ci)/3 + 0(e 2 ) (25) 



i.e. a — 1 — §e for a non periodic interface, and a = 1 



2 ti , „ , _ , 

for a charge density wave (CDW). The large-cutoff scale 
is given by v m with r\ m ~ m 2 ~ z , z = 2 — |e for non- 
periodic disorder and z = 2 - § for CDW [TMH] . 

To conclude, we introduced a general method to com- 
pute both spatial and temporal velocity correlations in 
an avalanche. Its tree-approximation is exact at and 
above the upper critical dimension d > d uc . There the 
center-of-mass motion is equivalent to the phenomenolog- 
ical ABBM model. This establishes the range of validity 
of the latter. For d < d uc corrections are calculated in a 
controled expansion in e = d uc — d. 
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